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TECHNICAL MEMORANDUM X-64886 

ON CONDITIONAL SAMPLING FOR TURBULENT FLOW STUDIES 

I. INTRODUCTION 


A conditional sampling technique has been used recently by many researchers in 
studying turbulent shear flows 1 1-9] . Under the conditional sampling, a discriminating 
criterion is chosen and flow measurements are made using two or more probes. One of the 
probes is the triggering probe and the rest are sampling probes. Averaging at the sampling 
probe is made only when the discriminating criterion is satisfied at the triggering probe. 

By so doing, the flow field is divided into different regions and average flow quantities in 
each region are obtained. A tremendous amount of data may be generated and by using 
the various forms of correlations, the detailed structure of the turbulent flow field is 
obtained. This is usually shielded from us in a conventional measurement. 

Any detectable flow quantity may be used to provide the discriminating criterion. 
The location of the turbulent-nonturbulent interface at the outer edge of a turbulent 
boundary layer has been used by Kovasznay et al. [ 1 ] , by Imaki [2] , and by Antonia [4] . 
The same has been used by Wagnanski and Fielder [5] in a two-dimensional mixing layer, 
and by Thomas [6] in a two-dimensional wake. The location of the viscous sublayer was 
used in studying turbulent “bursting” phenomena by Narahari et al. [7] and by Blackwelder 
and Kovasznay [8] . Other quantities that may be used include the flow temperature used 
by Dean [9] , the species concentration, and the location of a fluctuation shock front. 

Mathematically speaking, the conditional average is just a form of weighted average, 
with the weighting function generated from the reading at the triggering probe. By using 
different types of weighting functions, different types of conditional averages can be 
obtained. The conventional average, for which no discriminating criterion is used, becomes 
a special case. 

In the following sections, the concept of conditional sampling is demonstrated by 
considering the turbulent boundary layer over a flat plate. The flow field consists of two 
distinctive regions: the turbulent region near the flat plate, and the nonturbulent region of 
the ambient flow. These regions are separated by a thin, sharp interface that is highly 
contoured, and its position fluctuates in time. By using the location of the interface as a 
discriminating criterion, the following types of conditional averages can be obtained: the 
turbulent zone average, the nonturbulent zone average, and the point conditional average. 

A “second averaging” relation is then obtained which relates the conditional 
averages to the conventional average through the statistical distribution of the interface 
position. With this relation it is thus possible to incorporate the experimental data of 
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conditional measurement into the theoretical modeling of turbulent shear flow. Simplified 
relations for the “second averaging” are obtained under special assumptions. 


II. THE WEIGHTED TIME AVERAGES 


Consider the turbulent boundary layer over a flat plate as shown in Figure 1. 



Figure 1 . Turbulent boundary layer near the wall. 

Let the instantaneous location of the turbulent interface, assumed to be single- 
valued, be given by y = h(x,t). Let be a flow quantity measured at a point 

P(x,y). Then at any time t; one of the following conditions must be satisfied: 


h(x,t) > y , (la) 

h(x,t) < y , (lb) 

or 

h(x,t) * y , (lc) 

where the point P is said to be either inside the turbulent region, inside the nonturbulent 
region, or at the interface, depending on whether condition (la), (lb), or (lc) is satisfied. 

If a weighting function w( r ,t) can be defined that takes different values for P in 
different regions, then a weighted average of Q is defined as 
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K f~t\ - l»ni 0 

Q w (r) - 


/ Q(7, t) w(7,t) dt 


T 

J w(7,t) dt 
0 


According to the weighting function w( r ,t), the types of averages listed as follows 
can be obtained. 

1 . Conventional Averages. If a constant weighting function is used for P in all 
the three regions in (1), then the conventional average, or the “long time” average, is 
obtained, 


T 

/ Q(7,t)dt 

a?) • ± -j — 

/ dt 

o 


The fluctuation of Q ( r ,t) about its average Q ( r ) is given as 


q( r ,t) = Q( r ,t) - Q( r ) 


which has a “long time” average of zero, 


q(r) = 0 


The standard deviation of Q(r,t) is given by 


“ x “ Vi 

f ■ q(7 ,t) 1 dt 
lim 0 

o(r) = x-H>o ^ 

/« 

0 

which is a measure of the intensity of the fluctuation. 
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2. Zone Average. If the functional for the weighting function to be used is the 
turbulent signal 1 ( r ,t), which is unity in the turbulent zone and zero elsewhere, i.e. 



h(x,t) > y 


h(x,t) < y 


( 7 ) 


then the following turbulent zone average is obtained: 



lim 

T-*®o 


T 

/ Q(”r ,t) 1(7, t) dt 
0 


1 

/ 


I(r,t) dt 


( 8 ) 


A nonturbulent zone average can be obtained similarly by using 1-1 (r,t) as the weighting 
function, giving 



T 

/ Q(7,t) [1 - I(7,t)] dt 

lim __2 

T-+00 x 

f [l- 1(7, t)l dt 

Q 


( 9 ) 


Introducing an intermittency factor, 7 (r), as the conventional average of the 
turbulent signal, 
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/ I(7,t) dt 

*7) - 1(7) . *». JL_ , ( ,o, 
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then the turbulent and nonturbulent zone averages can be rewritten as 


/ Q(?,t) 1 ( 7 , t) dt 

• lim J> 

7 T ^°° T 

/ «« 

0 


/ Q(?,t) [1 - 1(7, t)] dt 


\ lim 0 
1 _ 7 T-k» 


T 

/* 

0 


Obviously, the conventional average can be obtained from the zone averages using the 
relation 


Q(r) = 7 (r)Q(r) + [1 - T(r)) Q(r) 


which was first obtained by Townsend [10]. 


3. Multiple Zone Averages. We shall now deduce from the zone averages the point 
conditional averages. First let us introduce the multiple zone averages. Assume that the 
range of the possible values of the interface position h(x,t) is given by the open interval 
(a, b), where a may be minus infinity and b may be plus infinity. Assuming that a 
sequence of ordered points Pq < Pj < Pj < ... < Pj exists such that « ‘a and 

Pj = b, then the interval (a,b) is divided into J smaller intervals by the sequence of 
points. If the following weighting function for each interval is defined 


1 , when pj_j < h(x,t) < pj 


Wj(r,t) = 


0 , otherwise 


for j * 1,2, 3, ..., J 
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then a multiple zone average at the jth interval is obtained as 


'•w — ► 

Qj(r) 


T 

/ Wj("r ,t) dt 

lim 0 
T-*» X 

/ w j (‘r,t)dt 

o 1 


t 


(15) 


Similar to equation (13), the conventional average is related to the multiple zone 
averages by the summation relation 



T 

Q(r) = £ fj(r ) Qj( r) 
j=l 


where 



_ lim 

- T-^oo 


T 

/ Wj(ir,t) dt 

0 ' 

T 

/ d. 

0 


(16) 


(17) 


is the probability that the interface position h(x,t) falls in the interval between pj_j 
and pj . A normalization condition has to be satisfied by these probability functions, 
which yields 
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£ fjtf) =1 08) 
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4. Point Conditional Averages. Let us now take the limit such that the largest 
interval in the multiple zone average o* item 3 becomes infinitesimally small, so that the 
weighting function in (IS) becomes the Dirac delta functions, 

1 , h = pj 

Wj(?,t) = 6(h - pj) = . (19) 

0 , otherwise 

Then the following point condition average is obtained: 



!_ 


I 


Qj(r) 


T 

/ CK? ,t) 5(h - pj) dt 

lim 2 

T-*» t 

/ 6(h- Pj )dt 
0 


( 20 ) 



In the same limit, as the largest interval goes to zero, a probability density function 
Pj(h) for the turbulent interface position can be introduced in the expression 

fj(r) = Pj (h) dh , (21) 


Substituting (21) into equation (16) and taking the same limit, the following } 

integral expression is obtained: s 

I 


00 

Q (t) = / Pj(h) Qj(r ) dh 

-OO 


(22a) i 


i 

t 

Here integration is taken over all possible values of h , assumed to be from minus f 

infinity to plus infinity. This can be rewritten as j 
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Q(r) * 


/ Pj(h) Qj(?) dh 

«oo 

oo 

/ Pj(h) dh 


(22b) 


by using the normalization relation 


/ Pj(h) dh = 1 


Noticing the analogy between equations (2) and (22), the relation between the con- 
ventional averages and the conditional averages can thus be interpreted as another weighted 

A ^ 

averaging process. Since the conditional average Q (r ) is already an averaged quantity, 
this relation is thus Mmed a “second averaging” relation. The weighting function for this 
“second averaging” is *ne probability density function P (h) of the turbulent interface 
position, which may take any value from 0 to . 


III. CORRELATIONS IN CONDITIONAL AVERAGES 

In this section, we shall derive the relations between the fluctuation relations 
in the conditional and conventional averages using the “second averagir Jon, equation 
( 22 ). 

Let the flow quantity to be averaged, Q(7,t), be the product o /o measurable 
fluctuating quantities, R( r ,t) and S(r,t);i.e., 

0(7,0 = R(r,t)S(?,t) (24) 


Here the dependency of Q, R, and S on the space cc dinate r is well under- 
stood and shall be omitted from the following discussions. 

Substituting equation (24) into equation (22) gives the relation. 


RS - J Pj(h)RSj dh 
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Separating R(t) and S(t) into their conventional mean and fluctuating parts gives 


R(t) S(t) = [R + r(D] [S + s(t)] ■ R S + R s(t) + S r(t) + r(t) s(t) 


Defining fluctuating parts in the point conditional average as 


gives similarly 


A A A 


R(t) S(t) = RS + Rsft) + Sr A (t) + r A (t)s.(t) 


Substituting (26) and (28) into (25) and using the relations 


r = s = r. 


= s A = ° 


RS + rs = / Pj(h) (R S + Oa>j 


which can be rewritten as 


FT 55 j Pj(h) (R S + r^A)j dh - / Pj(h) Rj dh f Pj(h) Sj dh 
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Thus the conventional correlation rs can be obtained from the point conditional 

A A 

averages R, S and the point conditional correlation r A s A once the probability density 

distribution function Pj(h) is known. Contracting equation (31) gives the fluctuating 
level, which is 

? - f Pj(h) (R 2 + £ ). dh - I / Pj(h) Rj dhl (32) 

-oo \-oo / 
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Higher order correlations can also be obtained in a similar manner. The expression 
always includes correlation of the same order in the ^oint conditional average and correla- 
tions of lower order. A third order correlation is given here as an example: 


- AAA A /S A A /\ 

qrs = / Pj(h) (Q R S + Q r A s A + R q A s A + S q A r A + q A r A s A ) j 


dh 


- QRS-Qrs-Rqs-Sqr 


(33) 


The above expressions can also be applied to two-point correlations like 


q(x,t) r(x + dx,t). Similar relations for time-delayed correlations like <q(x,t) r(x, t+dt)> 
can be obtained using ensemble averages. 


IV. SIMPLE EXAMPLES 
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Application of the “second averaging” relations, equations (22) and (3 1 ) through 
(33), in a real flow situation may become impractical when extensive conditional measure- 
ments are needed. Under certain conditions, however, approximations for the conditional 
averages can be made, and simplified relations for the second averaging can be obtained. 

In the following examples, similarity assumptions are used to draw up simplified relations. 


Example 1 . Fluctuations in an Infinite Domain. Consider the motion of a 
fluctuating interface in an infinite domain as shown in Figure 2. This interface may be a 
material surface which separates two different flow regions, or a nonmaterial surface such 
as a shock front or a shear front. It is reasonable to assume that the motion of the fluid 
very near the interface is always in phase with the motion of the interface; and the flow 
properties at a fixed distance from the interface position always remain constants, so that 
they can be given a single-valued function F as 
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which reduces to 


OO 

Q(x) * f P(x - y) F(y) dy 

-OO 


(35b) 


by using new variables 


y = x - h 


A particular example of this case is the averaging of the turbulent signal, which is 


F(y) = I(t) = 


0 

1 


y > 0 
y < 0 


(36) 


The average from (35b) gives 


0 °° 

I(x) = J P(x - y) d‘y - / P(h) dh = y 

-00 x 


(37) 


which gives the intermittency factor. 


Example 2. Linear Fluctuation Between Nodal Points. Consider two layers of 
mutually immiscible fluid bounded on both sides by solid walls as shown in Figure 3. 
Assume a one-dimensional fluctuating motion is produced in the x-direction, normal to 
the interface, and no reflective waves are produced so that the motion of the fluid is in 
phase with the motion of the interface. Furthermore, assume each fluid layer is stretched 
and compressed linearly by the motion of the interface. Then the point conditional 
averages in the fluid region can be represented by one of the two similar profiles given 
below: 



0 < x < h(t) 


h(t) < x < L 


(38) 
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The functions F, and F, are both smooth and single-valued, and they are normalized 
with the distance between the interface and the walls. 



Figure 3. Linear stretching of two fluid regions between two solid boundaries. 
Substituting equation (38) into equation (22) gives 


?; 


t 

it 


X 

Q(x) = / P(h)F 2 (j^.)dh + / P(h)F,^ 


dh 


(39) 


which, after changing variables, becomes 


Q(x) 




-I 


+ >, 


, F,(y) 


x 

L 


(y) dy ' 


(40) 
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The conditional sampling averages have been analyzed as weighted time averages. 
The various forms of condifioaai averages are obtained by using different types of 
weighting functions. A relation is derived which relates the conventional averages to the 
conditional Averages through a “second averaging” process. This “second averaging” 
process is a'so a weighted averaging over possible values of the discriminating quantity. 

A few examples are given in which simple expressions for the point conditional averages 
can be used. 

The “second averaging” relations, (22b) and (31) through (33), relate the condi- 
tional averages to the conventional averages. Application of relation (22b) in the mathe- 
matical modeling of turbulent shear flow has been done by Nee and Kovasznay [11] and 
by Wang [ 12] . Using a simplified form, encouraging improvements in the mean flow 
fields were obtained. The justification for these relations lies in the assumptions made on 
the conditional averages. Detailed conditional measurements are thus needed for this 
purpose. Recently, Paizis and Schwarz [13] used a rack of 20 probes in studying the 
structure of the turbulent interface. From the simultaneous measurements at the 20 
locations, detailed structure and motion of the interface were obtained. This information 
is valuable in the constructing of a better model for the turbulent shear flows. 
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